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Abstract. We present hybrid PIC simulations of fast electron transport and energy 
deposition in pre-compressed fusion targets, taking full account of collective magnetic 
effects and the hydrodynamic response of the background plasma. Results on actual 
ignition of an imploded fast ignition configuration are shown accounting for the 
increased beam divergence found in recent experiments [J.S. Green et al., Phys. Rev. 
Lett. 100, 015003 (2008)] and the reduction of the electron kinetic energy due to 
profile steepening predicted by advanced PIC simulations [B. Chrisman et al. Phys. 
Plasmas 15, 056309 (2008)]. Target ignition is studied as a function of injected electron 
energy, distance of cone-tip to dense core, initial divergence and kinetic energy of the 
relativistic electron beam. We found that beam coUimation reduces substantially the 
ignition energies of the cone-guided fuel configuration assumed here. 



1. Introduction 

Fast ignition [1], [2] involves transport of GA currents of laser-driven electrons through 
dense coronal plasma of imploded fusion targets. Recent studies |3l H] have shown 
that the beam can undergo the resistive filamentation instability when passing through 
coronal plasma and that beam collimation by self-generated fields can reduce ignition 
energies substantially. These results were obtained for rather focused beams (initial 
divergence half- angle = 22°) with electron kinetic energies given by the ponderomotive 
scaling. Recent experiments O |6] have evidenced an increase of beam divergence with 
laser intensity, e.g. electron effective propagation angles of 35° have been reported for 
fast ignition relevant conditions. In addition, advanced PIC simulations [7] have shown 
that the ponderomotive scaling overestimates fast electron kinetic energies due to the 
plasma electron density profile steepening caused by photon pressure at laser irradiances 
of the order of 10^° W/cm^. In this paper, we broaden the scope of our previous studies 
[3], H] and present results on actual ignition of an imploded fast ignition configuration 
accounting for the increase of the beam divergence and the reduction of the electron 
kinetic energy mentioned above. In addition, we assume a supergaussian distribution 
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in radius of the beam electrons to enhance azimuthal magnetic field generation and 
beam collimation. An imploded cone-guided fuel configuration is considered, and target 
ignition is studied as a function of injected electron energy, distance of cone-tip to dense 
core, initial divergence and mean kinetic energy of the relativistic electron beam. The 
present study has been motivated by the future fast ignition facilities such as HiPER 

In the hybrid approach pursued here, the relativistic electron beam is treated by 
2D/3D PIC including collisional energy loss, while the high density background plasma 
is modelled by resistive MHD equations including hydrodynamic motion to describe 
magnetic field suppression by plasma return currents. Full scale kinetic simulation may 
become possible in the future |7], but presently the hybrid approach pursued here offers 
a unique option to investigate important transport features such as current filamentation 
and magnetic beam collimation simultaneously with ignition physics (fusion reactions, 
a-particle transport and deposition, thermal radiation transport, hydrodynamics and 
heat conduction). One may recall that, so far, most fast ignition simulations [9] assumed 
ballistic straight-line beam transport, neglecting all the intricacies of high- current (GA) 
transport in plasma. 

The present paper is organised as follows. In Section 2, the simulation model and 
the initial distribution function assumed for the fast electron beam are described. In 
Section 3, the imploded fuel target configuration and the beam parameters used in our 
study are presented. In Section 4, we analyse the fast electron energy deposition in 
the target. In Sections 5 and 6, the coupling efficiencies and the ignition energies are 
parametrically studied, respectively, as a function of cone-tip to dense core distance, 
initial divergence and mean kinetic energy of the fast electron beam. Finally, in Section 
6 we summarise our results. 



2. Simulation model 



We used the hybrid model proposed by Bell [10] and further developed by Davies 
[TT] and Gremillet et al. [I2]. This model is adequate for describing self-magnetized 
transport in high-density fuel, where kinetic energy transfers and most of the beam- 
plasma instabilities are suppressed by collisions. It treats only the relativistic beam 
electrons by PIC and models the background plasma by the return current density j^, 
tied to the electric field E by Ohm's law with resistivity rj. Maxwell's equations are 
used in the form 

E = Vir, (1) 

jb+jr = — VxB, (2) 

dB/dt = - V X E, (3) 

where j=jf,-l-jr is the net current density, defined as the sum of beam and return current 
densities. The displacement current and charge separation effects can be neglected since 
in this high-density environment relaxation times and Debye lengths are much smaller 
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than the sub-picosecond and micrometer scales in which we are interested here. The 
recent study by Kemp et al. [13] justifies to choose resolution less than the Debye and 
skin lengths at high densities due to the suppression of collective effects by collisions. 
In the present model, the beam deposits energy into plasma electrons by direct classical 
Coulomb deposition and via return current Ohmic heating with power density r]j^. Our 
hybrid-PIC approach reproduces cone-target experiments p3l [15], similar to previous 
hybrid simulations [HI ITT] . Since laser pulses of about 15 ps are needed to ignite the pre- 
compressed fuel, we have included in our model thermal electron conduction, multigroup 
radiation transport, fusion reactions, alpha particle transport and energy deposition, 
and hydrodynamics. This integrated model is used here in 2D r - z cylindrical geometry 
mode. It is worth pointing out that previous hybrid simulations [H [161 HZ] assumed a 
stationary background plasma while in the present model plasma motion is described 
by radiation-hydrodynamic equations. The details of the fast electron transport code 
can be found in [HI [HI [20] . 

We assume in the present calculations that beam filamentation at the cone tip 
[T7[ [2T] or Weibel and two-stream instabilities [221 [23] in the low density zones do not 
perturb beam propagation from the critical surface to the minimum density of 2 g/ cm^ 
considered in this paper. This is supported by the cone target experiments carried out 
by Kodama et al. [HI [H] and also by the recent PIC calculations by Chrisman, Sentoku 
and Kemp [7] showing that, even though the fast electron beam is split into filaments 
at the cone tip, these filaments are strongly damped in the dense plasma beyond the 
cone. 

Because the ultra-intense laser pulse interaction with the cone cannot be described 
by hybrid codes, we rather model the injected beam in form of a directed energy 
distribution typically fitted to PIC simulations [3, |2l] or experiments [HI [HI [25] . The 
energy distribution of beam electrons is assumed to be one-dimensional (ID) relativistic 
Maxwelhan with temperature kTb = /t x mgC^ [(1 + /A^/1.37 x 10^^)^/^ — 1] obtained 
for local laser intensity I (in units of W/cm^) and wavelength A (in /im). Here, /t 
is a parameter that accounts for the plasma electron density profile steepening due 
to the laser ponderomotive force reported in Ref. [7]. We assume /t = 0.80, 0.67 
and 0.51 to obtain electron mean kinetic energies of 2.0, 1.6 and 1.2 MeV for a laser 
intensity of 2x10^° W/cm^ and A = 0.53 //m, respectiveljlil. We have chosen the ID 
relativistic Maxwellian distribution due to the ID character of the electron acceleration 
by ultra-intense laser pulses. In addition, the lack of experimental data about spectra 
of electron beams produced by laser pulses of a few times 10^° W/cm^ and 10 - 20 ps 
duration also supports our choice. It is worthwhile pointing out that the mean energy of 
the ID relativistic Maxwellian distribution {{E)i£, = kT in the ultra-relativistic limit) 
is substantially lower than that corresponding to the 3D relativistic Maxwelian one 
{{E)^!:) = 3kT in the ultra-relativistic limit) for the same laser intensity, leading to 
beam electrons with lower kinetic energy and penetration depth. 

J For electron temperatures kT;, = 2.24, 1.88 and 1.43 MeV, the ID relativistic Maxwellian distribution 
gives (£')/kTfc w 0.88, 0.86 and 0.84, respectively. 
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Figure 1. Energy deposition isocontours of (a) monoenergetic 2 MeV electrons 
impinging perpendicularly on a DT slab of 400 g/cm'^, and (b) electrons with the 
spectral distribution explained in the text, mean energy (E) — 2 MeV and divergence 
half-angle (9) = 22° (HWHM). Colour isocontours are in logarithmic scale in units of 
J/cm"^. (c) Comparison of the electron penetration depth obtained by different models. 

Initial divergence angle of the relativistic electron beam is a crucial parameter for 
fast ignition. Recent experiments of hot electron generation and transport in solid foil 
targets with IX^ (~ 4 x 10^^ W/im^/cm^ ) and pulse durations r (~ 5 ps) relevant 
for fast ignition have evidenced beam divergence angles of about 35° O [6]. Our 
hybrid simulations can reproduce this result by assuming an initial divergence half- 
angle (6) around 50°. The difference between this initial half-angle and the effective 
propagation full-angle of 35° is due to the beam coUimation by self-generated azimuthal 
magnetic fields Bg, which strength increase with pulse duration as shown by equation 
(3). Integrated simulations for initial half-angles (6) > 50° give electron-beam ignition 
energies Eig > 50 kJ, well over the energies envisaged for the future facilities [H]. Fast 
electron generation in cones may reduce the beam divergence due to surface acceleration 
and focusing toward the cone tip [26l [27]. However, there are experimental evidences 
showing almost no effect of the cone on beam focusing [28]. Here, we have considered 
the initial beam opening half-angle as a parameter, varying from the 22° measured 
by Kodama et al. [Ml |15] to the 50° used to reproduce the experiments mentioned 
above. The initial beam divergence is modelled as follows: a fast electron with kinetic 
energy E = {'j — l)meC^ is randomly injected in a cone with a half-angle given by the 
ponderomotive scaling formula 9 = tan^^j/g x [2/(7 — 1)]^/^}, where 7 is the Lorentz 
factor of the electron and the parameter fg is used to adjust the initial divergence half- 
angle of the whole electron beam. We have chosen fg = 1 for the reference case that 
corresponds to (6) = 35°, well below the initial divergence required to reproduce the 
experiments of Refs. [SlE], but substantially higher than the (6) = 22° considered in 
previous studies [3]- 

Fast electron collisions with the background DT plasma are modelled by the 
relativistic Fokker-Planck equation with standard Coulomb cross sections. The Fokker- 
Planck equation is solved by the Monte Carlo method [29]. The energy deposition of a 
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Figure 2. Central cut through imploded target configuration: (a) density isocontours 
in g/cm^ and (b) density profile at r = 0. 



monoenergetic 2 MeV electron beam with a supergaussian profile in radius and a half 
width at half maximum (HWHM) of 20 fim impinging perpendicularly on a 400 g/ cm^ 
DT slab is shown in figure [TJ The beam is coUimated at half of its range, being subject 
to scattering and blooming in the second half, as reported in [30]. The fast electron 
penetration depth obtained with the present model is similar to that recently published 
by Solodov and Betti [3T] and also to the quantum calculations by Starikov and Deutsch 
[321 [33] and about 20% larger than the scaling pR [g/cm2]= 0.6 xEhot [MeV] [31], as 
shown in figure [Tj^c). It is worth noticing that the energy deposition changes significantly 
when beams with realistic energy and angular distributions are considered. For instance, 
in figure [Tj^b) the energy deposition is shown for the same beam as in figure [Tj^a), but 
with a mean divergence half-angle of 22° (within the HWHM of the super-Gaussian 
profile in radius) and a ID relativistic Maxwellian distribution with mean energy (E) 
= 2 MeV. Notice the apparent range lengthening due to the beam energy spectrum and 
the subsequent delocalization of electron energy deposition. 

Regarding the energy transport models, we have used one-group flux-limited 
thermal electron conduction and multigroup radiation transport. Classical Spitzer 
electrical and thermal conductivities are chosen for the DT plasma with the Coulomb 
logarithms given in Ref. [35]. Degeneracy effects have not been taken into account in 
the transport coefficients because the DT is heated to keV temperatures in a few ps and 
field generation is not important in the dense core, as will be discussed in Section 4. We 
have neglected also the effects of self-generated magnetic fields on transport coefficients. 
MPQeos tables [36] have been used to compute electron and ion temperatures from 
the deposited energy. Alpha-particle transport and energy deposition calculations have 
been performed by a 3D Monte-Carlo model with standard stopping power and mean 
deflection coefficients. 
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3. Imploded target configuration and beam parameters. 

The target configuration used in tliis paper is sliown in figure O It is a simplified version 
of that obtained from radiation-hydrodynamics calculations |371l38]. It consists of 0.18 
mg of DT fuel compressed into a supergaussian spherical blob with a peak density of 
500 g/cm^ and 82 fim diameter (FWHM) placed on a density pedestal of 2 g/cm^. The 
pR of the blob is 2 g/cm^ and the pL through the central cut of the simulation box is 4 
g/cm^. For simplicity, a uniform initial DT temperature of 300 eV has been assumed, 
which sets the initial resistivity to a level of 3 x 10~^ Qm. We assume also that the DT 
is initially at rest, which is a reasonable choice because the fuel is almost stagnated at 
the time of peak pR. 

A beam of fast electrons is injected from the left at 2; = 0. We imagine that 
it emerges from the tip of a cone at this position with a supergaussian profile in 
radius and Gaussian in time. The supergaussian profile in radius reduces the beam 
energy spread and enhances azimuthal magnetic field growth because from equation (3) 
dBg/dt ~ dEz/dr ^ d{r]jb,z)/dr. However, the lack of experimental evidences of laser- 
generated electron beams with such profiles makes our calculations slightly optimistic in 
a similar way to other fast ignition studies that assume cylindrical beams [HI [M] . More 
realistic calculations with a Gaussian profile have been presented in Ref. 

The reference case considered in this paper is defined by a distance between the 
cone tip and the dense blob d = 100 pm, beam power of 2 PW, beam radius of 20 pm 
(HWHM), pulse duration of 18 ps at full width at half maximum (FWHM), electron 
mean kinetic energy within the HWHM of the radial distribution (E) = 1.6 MeV (/t 
= 0.67) and initial beam divergence half-angle (6) = 35° within the HWHM (/e = 1). 
The total beam energy is 36 kJ. Assuming a laser-to-fast electron conversion efficiency 
of 40%, this beam could be generated by a laser pulse with the same supergaussian 
distribution and spot radius that the electron beam, a mean intensity within the FWHM 
of 2.7x10^° W/cm^, A = 0.53 pm and a total energy of 90 kJ. The second harmonic 
of the neodymium laser is used here to reduce the mean energy of fast electrons and 
maximize their coupling with the dense core. The number density of the electron beam 
at the injection surface is around 2 x 10^^ cm~^, below the relativistic critical density 
for JA^ = 7.6 X 10^^ Wpm^/cm^, and a factor of ~ 5 higher than the standard critical 
density for A= 0.53 pm. 

Concerning the numerical parameters, we have chosen a cell width of 1 pm in each 
coordinate, a time step of 3 fs, and a total number of injected particles of 4 x 10^. Free 
boundaries have been used in all simulations. 

4. Fast electron energy deposition 

We present in this Section our results on the interaction and energy deposition of the 
fast electron beam with the target for the reference case defined above. 

In the low density halo, the beam deposits its energy mainly as Ohmic heating 
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Figure 3. Self-generated azimuthal magnetic field Bq (a) at the time of peak power 
and (b) at the end of the pulse. The fast electron pulse has a Gaussian profile in time 
with a total energy of 36 kJ and a duration of 18 ps (FWHM). Beam parameters are 
d = 100 /xm, {E) = 1.6 MeV and (6') = 35°. Dashed and dotted circles show the 
initial position of the 250 g/cm^ and 100 g/cm'^ isocontours, respectively. The solid 
line depicts the position of a perfectly coUimated beam of 20 /im radius. 



by return currents. Plasma electrons are heated up to temperatures Tg ~ 10 keV 
much higher than the ion temperature Tj in a few picoseconds. This temperature is 
almost constant while the electron beam is on because the energy transferred to plasma 
electrons is balanced approximately by thermal electron conduction in the direction 
perpendicular to the electron beam. The DT resistivity at 10 keV 4 x 10~^° fim) 
is high enough for a GA fast electron beam to generate an azimuthal magnetic field Bq 
of hundreds of Tesla in 10 - 15 ps [Bq rdE^jdr ^ rrjjb^z/rb = TVh/i'^rl), where h 
stands for total current, for beam radius and r for pulse duration). This magnetic 
field grows up slowly due to the low plasma resistivity and collimates the electron beam 
in the second half of the pulse. Azimuthal magnetic fields of a few hundreds of Tesla 
are sufficient to coUimate most of the beam electrons. For instance, the Larmor radius 
of 1 MeV electron in a i^e-field of 250 T is around 19 /im, similar to the 20 fim beam 
radius assumed here. If thermal conduction were not taken into account, temperatures 
in the halo would be much higher, the Sg-field would saturate at levels lower than the 
hundreds of Tesla reported here and there would be no significant beam collimation in 
the halo, as it was shown in previous studies [3]. 

At intermediate densities 10 < p < 100 g/cm^, electron temperatures are lower, 
resistivities higher and electron-ion energy exchange is fast enough to equilibrate electron 
and ion temperatures Ti ^ Te. Higher resistivities lead to enhanced field generation and 
collimation of the relativistic electron beam after a few ps, as depicted in figure [31 It 
is worthwhile noticing that the -B^-field peaks in this zone and how its oscillations 
are attenuated towards the beam axis due to the increasing plasma density. The 
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Figure 4. Energy density in deposited by the fast electron beam for the case 

of figure [31 Upper row [(a) and (b)]: Coulomb energy deposition and ion temperature 
at the end of the pulse in the full simulation with self-generated fields on. Lower row 
[(c) and (d)] Coulomb energy deposition and ion temperature at the end of the pulse 
in the simulation with self-generated fields artificially suppressed. Dashed circles show 
the initial position of the 250 g/cm^^ density isocontour. 

enhancement of the beam colhmation with time can be observed by comparing the 
location of the peak Bg-fie\d, which is near the beam edge, in Figs. [3](a) and (b) . 

At densities higher than 100 g/cm^, self-generated fields are damped by collisions of 
beam electrons with the background plasma and DT heating is almost exclusively due to 
Coulomb deposition. As has been pointed out recently [3], anomalous electron stopping 
due to resistivity effects plays no significant role at such high densities. The overall 
energy balance of the reference case shows that 87% of the beam energy is deposited 
in the target via Coulomb collisions, 12% is carried by electrons that pass through the 
target and escape, 1% is deposited as Ohmic heating by return currents and around 
54% is deposited at densities higher than half of the maximum (p > 250 g/cm^). Even 
though only 1% of the beam energy is deposited as Ohmic heating, self-generated fields 
turn out to contribute to the Coulomb energy deposition in the dense core indirectly via 
resistive colhmation of the fast electron beam, which concentrates the energy deposition 
in a spot slightly bigger than the initial beam diameter, as shown in figure ll](a). The 
effective propagation full-angle of the beam electrons near the dense core taken from 
Figs. [3] and m is, approximately, 28°, which gives a maximum spot diameter of 49 /xm 
(FWHM) at a depth of z = 72 fim. The energy density deposited via Coulomb collisions 
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Figure 5. Self-generated azimuthal magnetic field Be for the target pictured in figure 
[5] with d — 125 /im. The beam parameters are [E) = 1.6 MeV, [9) = 35° and a 
total energy of 50 kJ. Dashed and dotted circles show the initial position of the 250 
•gjcvc? and 100 g/cm^ isocountours, respectively. The solid line depicts the position of 
a perfectly coUimated beam of 20 [ivt\ radius. 



with and without self-generated fields are compared in Figs. |l](a) and (c). It is worth 
mentioning the enhancement of the energy density by a factor of about 2 when self- 
generated fields are taken into account. The oscillations shown in figure Hl^a) are due 
to the weak filamentation of the iJg-field in the density ramp. These oscillations do 
not affect DT ignition significantly because of thermal smoothing. However, for higher 
distances d from cone to dense core, the beam becomes more filamented (see figure [5]) 
and the energy deposition more fragmented than that shown in figure ID^a), imprinting 
the plasma density and temperature profiles. 

5. Coupling efficiency 

The electron beam coupling efficiencies (defined as the fraction of the beam energy 
deposited at densities higher than 250 g/cm^) predicted by our integrated model are 
depicted in figure [6t^a). The corresponding electron-beam ignition energies Eig are shown 
in figure int^b). In all simulations we keep constant the fast electron beam intensity and 
increase the beam energy firstly by increasing the pulse duration up to a maximum of 
20 ps and then by increasing the beam radius. For instance, the ignition energy of 57 
kJ is obtained for a pulse duration of 20 ps and a beam radius of 23.9 /xm while these 
parameters are 18 ps and 20 [xm for the reference case of 36 kJ. 

It is remarkable that, for the beam parameters chosen here, fast electrons propagate 
up to the top of the density ramp and deposit there a significant fraction of their 
energy without beam disruption or breaking-up due to the self-generated magnetic fields. 
Beams with lower kinetic energies are more prone to filamentation, while beams with 
higher kinetic energies have lower coupling efficiencies [3]. The importance of beam 
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Figure 6. (a) Coupling efHciencies and (b) minimum ignition energies of electron 
beams with (E) = 1.6 MeV as function of the cone-blob distance. Curves are labelled 
with the initial beam divergence half-angle (6) . The dashed lines and the lines labelled 
with 0° correspond to simulations with self-generated fields artificially suppressed. 

collimation can be seen in figure EJ^a) by comparing the coupling efficiencies for different 
d and {$) witli tliat obtained for a perfectly collimated beam {{$) = 0° and fields turned 
off). As shown in the figure, only for d = 75 fim and (9) = 30° the beam is nearly 
collimated, decreasing beam collimation with divergence angles, in agreement with the 
results of Bell and Kingham [39] . 

The dependence of the coupling efficiency on the mean kinetic energy of beam 
electrons is shown in figure [71 If self-generated fields are artificially suppressed, the 
optimal kinetic energy changes with the initial divergence half-angle from around 2 
MeV for perpendicular incidence to about 1.2 MeV for {6) = 35°. This is due to the 
different target areal densities seen by fast electrons with different divergence angles. In 
full simulations with fields on and (6) = 35°, the coupling efficiencies are much higher 
due to beam collimation and have a maximum around 1.2 MeV, similar to that found 
in simulations with fields off. It is shown also in this curve the weak dependence of 
the coupling efficiency on mean kinetic energy for (E) < 1.6 MeV. For higher kinetic 
energies, beam current density, self-generated fields and beam collimation are lower and 
the coupling efficiency tends to that found with fields off. 

Assuming a laser-to-fast electron conversion efficiency of 40% and a minimum 
overall coupling efficiency of the laser beam to the dense core of 25% |38], suitable 
fast ignition schemes have to have minimum electron beam coupling efficiencies around 
0.6. Figure El^a) shows that this coupling efficiency is found for {$) < 35° and distances 
d <100 /im, approximately, for the target analysed here. Figure [7] adds the additional 
restriction that (E) < 1.6 MeV to achieve overall coupling efficiencies around 25%. 
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Figure 7. Electron beam coupling efSciency as function of the mean kinetic energy of 
the beam electrons. Curves with solid markers correspond to simulations with Coulomb 
energy deposition only. The curve with open squares corresponds to full simulations 
including self- generated fields. Each curve is labelled with the initial beam divergence 
half-angle. In all cases the distance between the cone and the blob is d = 100 /im. 



6. Fuel ignition 

We have performed a series of integrated simulations to obtain the ignition energies 
Eig as a function of cone to blob distance d, initial divergence half-angle (6) and mean 
kinetic energy of fast electrons (E). Typical density and ion temperature profiles near 
the end of the pulse and in the ignition propagation phase are depicted in figure [H where 
one can see that DT expansion is not important during the beam energy deposition and 
becomes relevant after fusion reactions have been triggered. The ignition energies shown 
in this Section are obtained as the minimum beam energy for which the thermonuclear 
fusion power has an exponential or higher growth in time. Ignition energies as a function 
of d and (6) for a beam with (E) = 1.6 MeV are plotted in figure [6](b). This figure 
shows the very important dependence of Eig on both parameters and the crucial role 
played by self-generated fields. The ignition energies of figure Mjo) are 20-30% lower 
than those reported in previous studies |1] due to the sharper radial profile of the fast 
electron beam assumed here, i.e. supergaussian versus Gaussian, which concentrates the 
energy deposition and enhances magnetic field generation at the beam edge. Comparison 
of the ignition energies with that obtained for a perfectly coUimated beam {(6) = 0° 
and fields turned off) shows that for d = 75 /im and (6) = 30°, the beam is almost 
perfectly collimated. For higher d and (6), beam collimation is still important, but the 
beam diverges with an angle lower than the initial divergence assumed when propagates 
toward the blob. 

Sensitivity of the ignition energies on (6) and (E) is shown in figure [9] for a target 
with d = 100 fim. The relatively low sensitivity found for low initial divergences is a 
signal of beam collimation. For instance, the ignition energy Eig = 19.5 kJ obtained for 
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Figure 8. Ion temperature and plasma density isocontours for the case (a) and (b) 
of Figs. |4] with self-generated fields on (a) when 90% of the beam energy has been 
injected in the simulation box, and (b) 30 ps after the time of peak power of the fast 
electron beam when ignition is propagating through the DT fuel. Dashed and dotted 
circles show the initial position of the 250 g/cm'^ and 100 g/cm'^ density isocontours, 
respectively. Ion temperatures are given in keV and densities in units of 100 g/cm^. 



a beam with (E) = 1.6 MeV and (6) = 22° is lower than that obtained for a perfectly 
coUimated beam (22.5 kJ in figure [6t^b)) due to the beam compression by the i^g-field. 
It is also worth pointing out that the lowest ignition energies are obtained for (E) 
= 1.2 MeV, in agreement with the coupling efficiencies shown in figure [71 However, 
the dependence of Eig on (E) is more pronounced than that found for the coupling 
efficiency due to the variation of the penetration depth with the mean kinetic energy of 
beam electrons. 

Our integrated simulations predict that the imploded target configuration pictured 
in figure [2] heated by an electron beam with (E) = 1.6 MeV and {6) = 30° will ignite 
with beam energies E^g between 22 and 50 kJ, depending of the distance d between the 
cone and the blob. For the laser-to-fast-electron conversion efficiency of 0.4 assumed, 
these energies correspond to laser pulses from 55 to 125 kJ, respectively. Ignition of 
targets with d = 100 fim by laser beam energies lower than 100 kJ requires to generate 
electron beams with initial divergences and mean energies within the zone defined by 
Eig < 40 kJ in figure [9], which sets a compromise between the beam parameters (E) and 
(9) for the imploded target configuration analysed here. 



7. Conclusions 



One of the main conclusions of our study is that giga-ampere, multi-PW currents can be 
transported through the steep gradients of the plasma corona toward the high-density 
fuel core and can ignite it with beam energies of 30 - 40 kJ for the electron beam 




initial divergence half-angle (°) 

Figure 9. Electron beam ignition energies of the target depicted in figure [2] as a 
function of the initial divergence half-angle {0) with the mean kinetic energy {E) as a 
parameter. The distance between cone and blob is d = 100 //m in all curves. 

parameters and the target configuration assumed. The present simulations show that 
in the dense core, energy deposition takes place almost exclusively by classical Coulomb 
collisions. Self-generated fields play a major role for core heating improving the coupling 
efficiency substantially, but in an indirect way by means of beam collimation. 

We have performed a parametric study to obtain minimum ignition energies as a 
function of the cone - blob distance, beam divergence and electron kinetic energy for 
the imploded fuel configuration considered here. We found that the ignition energy 
depends strongly on the cone - blob distance and the beam divergence, both crucial 
parameters for fast ignition. The optimal kinetic energy for the initial divergence half- 
angle of 35° is about 1.2 MeV, significantly lower than the 2 MeV found for perfectly 
collimated beams. Our calculations show that targets with cone - blob distances lower 
than 125 /xm can be ignited by electron beams with energies around 40 kJ if initial 
divergence half-angles are about 30° - 35° and electron kinetic energies are lower than 
1.6 MeV. Assuming a laser-to-fast electron conversion efficiency of 40%, these electron 
beams could be generated by the short-pulse laser beams around 100 kJ envisioned for 
future facilities [8]. These conclusions rely on the integrated hybrid PIC model used 
here, which is valid for electron transport in dense media and therefore depends on 
the distribution function assumed for the injected electrons. Further theoretical and 
experimental investigations for a full characterization of fast electron beams generated 
in cones by laser pulses of 10 - 20 picoseconds are necessary to estimate more precisely 
the energy requirements of electron-driven fast ignition. 
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